Interaction of a dipole point vortex with flat boundary 
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Abstract 

In this work we have found an exact solution for the problem of the movement of a 
dipole type point vortex in an area of fluid limited by a flat boundary. We also present 
a solution to the problem of dipole point vortex motion in a right angle. It is shown 
that unlike a usual point vortex, the dipole vortex always comes away from the boundary 
asymptotically. This important feature of the dipole vortex allows it to be considered to 
be one of the efficient mechanisms of vorticity transfer from boundary to media 

1 Introduction 

The new type of dipole point vortex was found in work [1]. Also in this work, the 
motion equation system was obtained and studied for one-dimensional singularities which 
are compatible with the Euler equation. We point out that the point dipole vortices 
are weak solutions of Euler equation. It should be noted that point dipole vortices are 
Euler equation solutions only in cases where dipole moments themselves are evolving over 
time. This kind of vortex can be considered as a specific hydrodynamic quasiparticle as well 
as an already-recognized point vortex. The interaction between dipole vortices and usual 
point vortices defines a field of velocities in two-dimensional hydrodynamics. (Singularities 
of higher orders are not dynamically compatible with the Euler equation, but they can 
exist as stationary solutions.) The system of any number of such vortices is hamiltonian 
and has three integrals of motion in involution. According to the Liouville theorem, this 
means that the problem of interaction of a usual vortex and the other, of the dipole type 
[Ij is exactly integrated in a similar manner, like the known case of integration of system of 
three usual point vortices [2], [3], [4]. Exact solutions for, and the behaviour of two-point 
vortices, one of which is the point dipole, were considered in work [5]. In work [6] new 
stationary solutions with complex singuliarities of two-dimensional ideal hydrodynamics 
were obtained with help of point dipoles. But in these works, point dipole vortices were 
considered in fiuid without boundaries. In this work, the simplest case of the motion of one 
point dipole vortex with a fiat boundary to media is discussecl^. This problem is exactly 
integrable using the method of images and gives a complete description of all the modes 

^When our work was sent to the Journal we learned that two weeks earlier, work [9] had appeared in which 
motion equations of the dipole vortex near the plane boundary and solutions for dipole vortex movement were 
obtained, independently and using a different method (complex variables) In our work we have used the images 
method and analyzed in more detail the behavior of a dipole vortex near a plane boundary. Additionally, we 
studied the behavior of a vortex near a boundary which has the form of a right angle. This strengthens the 
main conclusion of our work concerning the role of a dipole vortex in efficient vorticity exchange between a 
boundary and media. 
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Figure 1: Position of the boundary and the point dipole vortex is shown. The fluid occupies 
semispace Xi > 0. 



of motion of a vortex near a flat boundary. In addition, we examine point dipole vortex 
motion at a right angle. As a result, we have established the simple asymptotic laws of 
vortex evolution on a long time scale and have found that the dipole vortex always moves 
away from the boundary over time. This means that we obtained an effective carrier of 
vorticity from boundary to media. This property is important, since, as a rule, vorticity 
origination is connected to boundaries. 



2 Search for the image 

To determine the movement of a dipole point vortex near a solid boundary it is possible 
to use the method of images [7], [8]. According to this method, the point dipole vortex at 
a wall corresponds to the system of two point vortices of a dipole type in fluid without 
a boundary. Such a system of point vortices must satisfy the same boundary conditions 
as the initial problem. Let the dipole vortex be at flat boundary with fluid xi = which 
fills the semispace xi > (see Fig{T]). With this geometry, the stream function must 
satisfy the boundary condition (p\xi=o = 0. Hence it is necessary to find the conditions 
when the dipole vortex, interacting with its image, would satisfy the boundary condition 
described above. The vortex-image settles at points x^^^ = —xi and Xg^^ = X2 for reasons 
of symmetry. It remains to find out how the dipole moment of the image D^^'f is connected 
with the dipole moment of the initial vortex D = [Di,D2). We can note the necessary 
conditions using the stream function of these two vortices: 

= (1) 

xi=0 

Here x{t) is the position of the initial vortex, and x{t)^^^ is the position of its image. Taking 
into account the relation between the vortices' coordinates, this condition transforms easily 
into 

-Di{t)xi{t) + D2{t){x2 - X2{t)) + D^^\t)xi{t) + D^^\t){x2 - X2{t)) = 
and it is satisfied when 

I?f)(t) = Di(t), D^^\t) = -D2it) (2) 

Then, the behaviour of a vortex near the solid boundary can be described in the boundless 
media with two point dipole vortices with the special choice of their characteristics. 

The component of the image dipole moment which is normal to the boundary coincides 
with the normal component of the initial vortex and the tangential image component 
reverses the sign. Using this fact one can consider the behaviour of the dipole vortex for 
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Figure 2: Shows the position of the dipole vortex in the point x = {xi,X2) with the dipole 
moment D = {Di,D2) in a right angle. The same boundary conditions appear when three 
other dipole vortices images are located in an unbounded media with the coordinates x^^^ = 
{—xi,X2), x^"^^ = (xi,— X2), x^^^ = {—Xi,—X2) and dipole moments chosen in accordance with 
= {Di,-D2), = {-Di,D2), = (^-D^,-D2). (The superscript in parentheses 
corresponds to the vortex number). 



other boundary conditions as well. In particular, it is easy to study the dipole vortex 
motion in a fluid limited by solid boundary forming a right angle ( the area occupied by 
fluid xi > and X2 > 0). 

The current function vanishing at the media boundary can be obtained by placing 
image vortices as shown on Figi2j In the case of the right angle boundary three image 
vortices are sufficient. It is easy to verify that the current function of this dipole vortices 
configuration vanishes with xi = and X2 = 0. Consequently, the dipole vortex evolution 
in a right angle comes to the problem of motion of four dipole vortices in boundless media. 
The values of these dipole vortices moments are shown on Figj2j In much the same way, 
it is easy to obtain vortices motion equations for angles of a different value as well as for 
a circular area and other elementary boundaries. 



3 Equations of motion 



To obtain the hamiltonian of this system of vortices we will use the general hamiltonian 
which was found in work [T]. In accordance with this work, the hamiltonian of two 
interacting dipole vortices is 

27r |x-x(i)|4 

Here x, x^^^ - are the coordinates and D, D^^^ - are the dipole moments of two dipole 
vortices. Now let us consider the relation which was demonstrated earlier between the 
positions and dipole moments of the vortex and its image. After simple transformations 
we obtain 

2tt 4xf ^ ^ 

Based on this hamiltonian it is easy to get the equations of motion of the dipole vortex 
near a solid wall in the form: 

dxi _ dH 



— = -e-k— (6) 
dt * dxk 

Where £ik is the antisymmetric unit tensor. After a substitution of the Hamiltonian 
we obtain the following system of equations: 



dxi D2 



dt AttxI 
dx2 Di 



dt Airxj 



(7) 
(8) 



dD^ 

(9) 



dt 

dD2 Dl + Dl 



(10) 



dt 47rx^ 

Firstly, from this system of equations follows the conservation of the orthogonal to the 
boundary component of the dipole moment of the vortex Di{t) = Di{0) = const. This 

is the consequence of the law of general conservation / = Yla^<^^vi ~ J2f5-^i^\ which 
is satisfied for the system of interacting usual point vortices and point dipole vortices 
[T]. In this equation is the vortex strength of a-point vortex and a;[,°''-its coordinates. 
Therefore the system splits into two subsystems of equations. It is enough to solve the 
closed system of two equations ([7]), to get all the vortex characteristics. It is important 
to note that the system of equations (|7l)- (ll0p can be directly obtained from the equations 
of motion of two dipole vortices [1] , after substitution of the relation between coordinates 
and the dipole moment^. It is evident that the energy of a point vortex is conserved and 
therefore the value of the dipole moment is defined by the distance to the boundary 
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where Eq is the initial value of vortex energy. Then it obviously follows that 

D\ + D\ = -SttEqxI 

We can define the dependence on time of the longitudinal movement to the boundary 
component of the dipole moment D2 using this invariant. In order to do it we solve this 
equation with respect to xi 

Here we take into account that energy of the dipole vortex E'o < and fluid occupies the 
positive semiplane x\ > 0. This defines the choice of the positive sign of the square root. 
We exclude coordinate xi from the equation (|10|) with help of this formula. As a result 
we get 

dD2 C 



dt 



where the constant C = (— SttE'o) ^ /47r > is defined by vortex energy. Taking into 
account the conservation of Di , we can integrate this equation in elementary functions 



D2\ Dl + Dl + Di\n[D2 + J Df + = 2Ct + const (12) 



^For the equation system ([7|- pIH) given above we used different variables than those in work [5]. 



Figure 3: Change of D2 in time: on the left for initial conditions -D2(0) = 1, -Di(O) = 1, 
xi(0) = 1 and X2(0) = 0, on the right 1^2(0) = — 1 and other initial data coincide with the 
previous initial conditions. 



The integration constant is defined by the initial conditions in accordance with the equa- 
tion const = 1)2(0)^^1+^2(0) + -^1 In (-^2(0) + V^i +^2(0)) • Here, and further in 
this article, for initial values of D2{t) we use the designation of the 1)2(0) type. From this 
solution it is easy to find the asymptotic increase of D2{t) according to the relation 

D2{t) ~ y/2Ci 

Naturally, at the given law of evolution D2{t), the change of coordinates of the vortex 
can be easily calculated. Then xi{t) is defined from D2{t) by an algebraic equation (lllh 
which can be written as 

,^? + ^i(o), 

using the law of energy conservation. After that, the coordinate is obtained by the in- 
tegration of the equation ([8]). The degeneration solution with Di = was considered in 
article [9]. In this case in accordance with equation (jl3p the position of the dipole vortex 
xi = D2 ^2(0) proportional to the non-zero dipole moment component, and the solution 
of equation ([JO]) takes the evident form D'^it) = D2{0)^ ■ (1 + i^^^^)- It follows that 
with 1)2(0) < the dipole momentum vanishes over a finite time and consequently only 
in this case the dipole vortex approaches the boundary. With any small deviation from 
zero of the normal to the boundary dipole moment Di component, the approach stage to 
the boundary switches after several times to move away from boundary to infinity. 



4 Movement of vortex at solid boundary 

Let us now consider the character of the movement of the vortex at the wall with help of 
the previously-obtained exact solutions. First of all, one can notice, that the component 
of the dipole moment along the wall increases with time. Over long times it follows the 
square root law. At the initial stage at 1)2(0) < this component decreases at first, and 
then increases (see FigJS]). 

As a result of this behaviour of the longitudinal component of the dipole moment, over 
long times the vortex moves asympoticaly away from boundary as 



2Ct 
D'iTDliO), 

Hence, for 1)2(0) < at the beginning, the vortex is approaching the boundary at a 
minimal distance (see eq. (fT3]) ) 
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Figure 4: On the left is shown the trajectory of vortex movement for Di > and D2 > and 
on the right for Di > and D2 < 0. The direction of movement is marked by the arrow on the 
trajectory. On the right, one can see the beginning stage of the approach to the boundary and 
further moving away from it. 



Then the vortex moves away from the boundary asymptotically in accordance with formula 
(I14p . Now we have to discuss the movement of vortex along the boundary. For this purpose 
we will come back to equation ([8]). The equations ^ and pT]) give: 

dxi ^ ^/{D'i + D'i{0))xi - xi{0)^Dl 
dx2 -DiXi(O) 

Integrating this equation we obtain the vortex movement trajectory along the boundary: 



X2 

-sign{D2{0))— j— + const 



In ^^D-i + D-i{0)x, + ^{D-i + Dj{0))x'i-D-ixi{0)' 

(15) 

The integration constant is defined by initial conditions. Examples of the trajectories of 
vortices' movement are shown in Fig. HI 

The main conclusion that follows from the equation (jlSp is that the vortex moves 
away exponentially from the wall. It is important to note that a usual point vortex moves 
with constant speed along a flat boundary maintaning the same distance from it. The 
direction of the longitudinal motion of dipole vortex is defined by the sign of -Di(O). From 
a physical point of view, this means that vortices of a dipole type can generate an effective 
mechanism of vorticity transfer from the boundary where it is generated to the media. 
Naturally, the dipole vortices can also intensify the transfer of other passive "impurities", 
such as temperature to media. For many physical phenomena this is a highly important 
property. 

In conclusion, we must highlight, that the case of two interacting dipole vortices is 
integrated in quadratures with the special choice of dipole moments. This choice is con- 
nected with zero values of some first integrals of motion. The case considered above 
corresponds to zero value of / = ^ D2 = and J = Yl • = 0. The remaining 
integral ^ Df = const and energy can take any value. In general cases, the problem of 
dynamics of two point dipole vortices is not integrated. 



5 Movement of dipole vortex in a right angle 

Using the images method one can obtain the hamiltonian of the dipole vortex in media 
limited by a right angle. With help of the hamiltonian of four interacting dipole vortices 
and after the substitution of coordinates and dipole moments of images we obtain: 

I \Dl + Dl Dl + Dl {xl-xl){Dl- Dl) + ADiD2XiX2 
n = — — < n 1- 
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Figure 5: Movement of the dipole vortex with the initial coordinates xi(0) = 0.8, X2(0) = 0.8 
and initial value of the dipole moment -Di(O) = 1.1 D2{0) = —1, which are obtained numericaly. 
A similar qualitative modification of vortex movement appears under any minor breaking of 
symmetry, for instance with conservation of Di = —D2, but with deviation from the diagonal 
xi(0) 7^x2(0). 

Here xi > and X2 > 0. As earlier, this hamiltonian can be considered to be like the 
hamiltonian of dipole vortex moving in an angle. We obtain the motion equation using 
the standard way in accordance with the equations: 

dxi _ J_ f £2 £2 D2{xl - xl) - 2DiXiX2 



dt 2tt \ x\ X2 {x\ + xlY 



dx2 _ 1 (Di^Di Di{xl - xl) + 2D2X1X2 



dt 2tt \ xf xl {xf + xl)'^ 

dDi 1 [Dl + Dl {Dl- Dl)x2 + 2DiD2Xi ^ (Dj - Dl){xl - xl) + AD1D2X1X2 
~ H r-T— — — 2x2 



dt 27r \ xl {xj + xl)'^ {xj + xl)^ 

dD2 _ 1 ( Dl + Dl {Dj - Dl)xi + 2D1D2X2 ^ {Df - Dl){x\ - xl) + AD1D2X1X2 

dt ~ 27t\ xl + (xf + Xi)2 {xl + x2)3 

Obviously this is a more complex dynamical system with two degrees of freedom. The 
energy of this dipole vortices system is conserved. However, other conservation laws 
obtained in work [T] become trivial. For the considered configuration of the dipole vortices 
they vanish. That is why the question regarding the integrability in quadrature of this 
equation system remains open. But it is easy to see the discrete symmetry of equations 
which is due to the permutation of coordinates. This permits us to find the simple 
particular case of dipole vortex movement with xi = X2 and Di = —D2. The motion 
equations for it are simplified and take the form: 

dxi 5 Di 
dt Att xl 

dDi _ 5 Dl 
dt ~ A-K x\ 

This particular case is exactly integrable. And the solution of the equation system has 
the form: 

'^i(0)~ 



Figure 6: Shows the trajectory projections of the movement of vortices on plane (xi, X2)with 
different initial dipole moments, as well as different initial positions of the dipole vortices in 
the media. On the left Di(0) = 1, ^2(0) = -1, and on the right Di{0) = -1, ^2(0) = l.On 
the left only one trajectory in accordance with the analytical solution obtained above, moves 
to the angle vertex. On the right, the same trajectory moves away from the angle vertex along 
the bisector. 



"^i-^^i^"^ 2^i(0) 

Hence, with -Di(O) > the dipole vortex approaches the angle vertex, and with initial 
condition -Di(O) < 0, moves away along the angle bisector from its vertex towards media. 
The character of dipole vortex movement is sensitive to deviations of the condition Di = 
— Z?2- So, even with small deviations of this condition, the dipole vortex approaching the 
angle vertex begins to move away from it over time (see FigH]) We can also examine the 
asymptotical behaviour of the vortex for example with X2 ^ xi. In zero approximation, 
the vortices motion equations system coincides with equations (f71)- (fT0l) . This allows us to 
describe the asymptotical behaviour of the vortex. The numerical results also demonstrate 
that the moving away of dipole vortices from angle vertex to media is typical (see. FiglGj). 
Consequently, for the case of a right angle, the dipole performs an efficient exchange of 
vorticity between boundary and media. 



6 Conclusion 

With the dipole point vortex the set of specific quasi-particle increases. This allows us 
to consider problems of bi-dimensional hydrodynamics within the framework of finite- 
dimensional dynamics systems. Hence the number of problems which are needed to take 
into account the impact of point dipole vortices increases as well. They rank with usual 
point vortices, but they also have some interesting properties. 

On one hand, dipole vortices are an example of unusual hamiltonian systems which 
appear without using Lagrange's formalism and Legendre's transformations. The hamil- 
tonian nature of the dynamical equations of interacting dipole vortices plays an important 
part in the search of equation solutions as well as for many other problems. New exam- 
ples of exactly integrable problems appear due to the existence of additional conservation 
laws. The examples of these exactly integrable systems are given above, but some others 
are known, for which references were given in the Introduction. On the other hand, these 
point dipole vortices have different properties compared to the usual point vortices. In 
particular, as was shown earlier with the solution of the problem of point dipole movement 
for an area limited by a plane boundary as well as for an area limited by a right angle, 
the important function of the point dipole vortices may be to provide a vorticity exchange 



mechanism between boundaries and media. This property is very important for many 
phenomena concerned with vortex movements, because as a rule, vorticity is generated 
near the boundaries, as in the classical examples of the flow around a body, appearance 
of von Karman vortex street, wing lift mechanism etc. From this point of view, it is 
interesting to consider the advection which appears in the field of interacting usual and 
point dipole vortices. It is evident also that there are more problems concerned with the 
statistical description of point vortices system. The statistical description of the singular- 
ities set stands in close relation to bi-dimensional turbulence. The appearance of dipole 
point vortices in superfluid liquid and in other complex hydrodynamical mediums is also 
of great interest, as well as the existence of point dipole vortices on the bi-dimensional 
manifolds. Most of these problems remain unanswered. 
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